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Abstract 

This paper investigates Automatic Generation Control (AGC) 
of hydrothermal system under deregulated environment. A 
novel dual mode controller is used in AGC and the 
optimization process is done by Fuzzy Adaptive Bacterial 
Foraging Method (FABFM). The dynamic performance of 
the system with regard to peak time, overshoot and settling 
time has been improved by means of Flexible AC 
Transmission Systems (FACTS) devices like Static 
Synchronous Series Compensator (SSSC) and Thyristor 
Controlled Phase Shifter (TCPS). There might be problems 
with present controllers to be solved by soft computing 
techniques. Some factors influencing the formulation of 
AGC problem are open transmission access and the evolving 
of more socialized companies for generation, transmission 
and distribution. Therefore, the effect of bilateral contracts 
on the dynamics should be considered in modern AGC 
systems. Results reveal that using dual mode controller 
based on FABFM plus SSSC and TCPS has much more 
dynamic efficiency than the system with integral controller. 
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Introduction 

The objective of a power system is to keep different 
kinds of balances such as between load-generation, 
scheduled and actual tie line flows. For the purpose of 
constant frequency as the primary index of healthy 
operation of system as well as the quality of supplied 
power to consumer, these two factors are employed. 
Considering load demand, generation is the parameter 
that maintains both of these. As frequency becomes 
low, generation is increased and if the actual outflow 
is greater than the scheduled, then it will be decreased. 
Since system conditions are always altering as a result 
of load variations during different hours of a day, it 
would be impractical to control these balances. The 
aim of AGC is to keep frequency almost constant and 
to regulate tie line flows (C. Concordia, 1954- M.L. 



Kothari, 1980- J. Nanda, 1983). The power system 
structure is changed, consequently, the growth of 
more specialized industries for generation (Genco), 
transmission (Transco) and distribution (Disco), under 
open market system (deregulation) has become 
conceivable. (R.P. Schulte, 1996) described the control 
of generation in deregulated power systems. To 
equilibrate reliability with economics Independent 
System Operator (ISO) as an unbiased coordinator has 
also been studied in (J. Kumar, 1997). The estimation 
of Automatic Generation control in a deregulated 
environment has been perused in [B. H. Bakken, 1998]. 
In addition, a detailed review of the subject and the 
simulation of an AGC system after deregulation is 
presented. Furthermore, making use of power 
electronic devices in power system control as FACTS 
results in flexible power system operation and control 
(P. Bhatt, 2010). In order to analyze the interconnected 
power system, SSSC and TCPS are considered 
beneficent for the tie-line power flow control of an 
interconnected power system. 

Lately, different intelligent technique have been used 
widely in different fields of electrical engineering. A 
fuzzy controller in a three -area system in the presence 
of SMES units on each area has been presented in (A. 
Demiroren and E. Yesil, 2004). In (G. Panda, S. Panda 
and C. Ardil, 2009), by means of fuzzy controllers, a 
two-area thermal system has been simulated on the 
existence of limitations in the generation rate of units 
and the results obtained from this method have been 
compared with that by common proportional integral 
controllers. In (P. Bhatt, R. Roy and S.P. Ghoshal, 2009), 
a two-area system has been studied in the presence of 
TCPS and systems of capacitor energy storage where 
the designed PID-controller parameters have 
optimally been regulated using various methods of 
artificial intelligence. In (R.J. Abraham, D.Das and A. 
Patra, 2006) the effects of applying TCPS on 
controlling the real power exchange between areas of 
generation control system have been examined. 
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This paper presents the design of a dual mode 
controller for AGC of hydrothermal system via 
FABFM. The dynamic performance of the system with 
regard to peak time, overshoot and settling time has 
been improved by means of FACTS devices like SSSC 
and TCPS. Open transmission access and the evolving 
of more socialized companies for generation, 
transmission and distribution affect the formulation of 
AGC problem. Hence, the traditional AGC system is 
modified to consider the impact of bilateral contracts 
on the dynamics. Simulation results have 
demonstrated better dynamic performance via the 
proposed controller compared with integral controller 
and without SSSC and TCPS. 

Problem Statement 

Dynamic Model 

In an AGC system, there is a connection between two 
areas under deregulated scenario. Area 1 consists of a 
reheat system and area 2 comprises hydro system. 
Here, in each area, three generators are supposed. 




FIG. 1 TWO AREA HYDROTHERMAL SYSTEM 
Fig. 1 depicts the block diagram of two-area 



hydrothermal system under open market scenario in 
which ACE of every area is fed to the corresponding 
controller. For every incoming ACE at the specific load 
variation, the precise control signal is produced. To 
compare the efficiency of both systems, a performance 
index is used as follows: 







(1) 



Modeling of SSSC 



In SSSC self-commutated voltage-source switching 
converters combines a three-phase voltage in 
quadrature with the line current. Such converters 
emulate an inductive or a capacitive reactance 
affecting the power flow in the transmission lines. The 
magnitude and polarity of injected voltage, Vs are the 
two factors for dynamically controlling compensation 
level. The device is in dual mode of capacitive and 
inductive. In order to stabilize the area frequency 
oscillations, the schematic of an SSSC in series with the 
tie-line between the interconnected areas by means of 
high speed controlling of the power is exerted, as 
shown in Fig. 2. The equivalent circuit of this system 
can also be represented by voltage source Vs in series 
with a transformer leakage reactance Xs. The voltage 
source Vs is the controllable factor of SSSC that 
actually represents the magnitude of injected voltage. 
Fig. 3 demonstrates the phasor diagram of the system 
considering the operating conditions of SSSC. 

Vm<em 




Area 1 



\ 



Area 2 



FIG. 2 SCHEMATIC OF SSSC APPLIED TO AGC 

Based on the above figure when Vs =0 , the current Io 
of the system can be written as: 



V -V 
j __ m n 

o 



(2) 



Where Xt =Xl+Xs. The phase angle of the current can 
be expressed as: 



=tan 1 

c 



-V cos0 
m m 



V sin& -V sinff 
m m n n 



(3) 
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FIG. 3 PHASOR DIAGRAM AT VS=0 AND VS LAGGING I BY 90 
But Eqn (1) can be expressed in a generalized form as: 

-V 



V -V -V 
m s n 



J X 1 



V -V 
m n 



J X 1 



I +AI 

o 



(4) 



The term AI is an additional current term due to SSSC 
voltage Vs . The power flow from bus m to bus n can 
be written as Smn=VmI*=Smno +ASm«, which implies 



1 inn 



+ jQmn ~ (Pmno + ^mn ) + J(Qmno + ^Qmn ) (p) 



Where Pmm and Qmno are the real and reactive power 
flow, respectively when Vs=0. The change in real 
power flow caused by SSSC voltage is given by: 



AP, 



inn 



Y^ s in(0 m -a) 

A. 



(6) 



When Vs lags the current by 90°, APm„ can be written 
as: 



AP, 



mn 



Or) 



(7) 



From Eqn (2) the term cos(9m-9c) can be written as: 



cos(0 m -0 C ) = -^-cos(0 n - C ) 



v, 



m 



Referring to Fig. 3, it can be written as: 



cos(0 n -0 c ) = 



yw 
xy 



(8) 



(9) 



And it can be seen as yw=Vmsin 0mn 



xy: 



I 2 2 

: V y OT + V n ~ ZVmVn cos0 mn and mn = m -0 n (10) 
Based on these relationships, Eqn (6) can be modified 
as follows: 



&P„ 



^Lsin(0, 



)> 



T ijv£+v£ -2V m V n cos 0„ 



(11) 



From Eqn (4), it can be written clS Pmn — Pmno +AP 
which implies 



T yjvA+Vii-2V m V n co S 9 mn 

Linearizing Eqn (12) about an operating point, it can 
be written as: 



%^cos(# m -0 n )(A0 m -A0 n ) 



+ (%^ S in(0 m)J )x- 



(13) 



A T JvZ+vZ -2V m V n cos 0„ 

l\ 1 mn— APtie+ APsssc which implies 



AP, 



SSSC 



%^sin(0 m „)x- 



AV, 



(14) 



L P V y m +V n -2V m V n cos0 mn 

As stated in Eqn (14), it is clear that variations in the 
SSSC voltage AVs, result in the controlling of SSSC's 
power which will consequently control the frequency 
and tie line deviations. The structure of SSSC to be 
included in the two area system in favor of reducing 
the frequency deviations is provided in Fig. 4. The 
frequency deviation of area 1 would act as input of the 
SSSC device. 



Afi(s) 


(1+sTi) 




(l-hiTj) 






(1+sTz) 

















1 


APsssc 




(1+sTsssc) 





FIG. 4 SSSC STRUCTURE 
Modeling of TCPS 

Fig. 5 shows a two-area hydrothermal interconnected 
system assuming TCPS in series with the tie-line. 




FIG. 5 SCHEMATIC OF TCPS IN SERIES WITH TIE LINE 

Area 1 is the thermal one comprising three reheat 
units is followed by TCPS. Second area is the hydro 
one consisting of three hydro units. The incremental 
tie-line power flow from first area to second one with 
TCPS under open market system is expressed as Eqn 
(16). 

The phase shifter angle A cf) (s) can be written as: 

-AError^(s) 



A<P(S): 



1 + sT 



(15) 



ps 



Where K4 and T ps are the gain and time constants of 
the TCPS and AErron(s) is the control signal which 
controls the phase angle of the phase shifter. Therefore, 
it can be written as: 



AP. 



tie\2 



(j) = 



^-[af 1 (5)-AF 2 (5)]h 



K n 



(16) 



Pn x T- AError, (s) 

12 l + rtps 1 
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AErron can be any signal such as the thermal area 
frequency deviation Afi or hydro area frequency 
deviation Af2 . 

Designing of Dual Mode Controller 

The steady state zero error could be achieved by just 
implementing the well designed integral controller, 
but the response of the system becomes slow resulting 
from high overshoot and long settling time. The 
proportional controller by itself will improve 
overshoot and the speed of the response. Apparently, 
the presence of the proportional controller is 
compelling at transient states to make system response 
faster so as to reduce the overshoot. The main 
shortcoming of the proportional controller is that it 
fails to bring the steady state error to zero. So the 
usage of both the integral and proportional controller 
simultaneously seems urgent. In this study, effort has 
been made to develop a dual mode controller 

Discontinuous Mode 

The control law employed during the transient period, 
i.e., the discontinuous mode, is switched between Eqn 
(17) and Eqn (18) depending on the magnitude of error 
signal i.e., ACE(t). For IACE(t)l>e, the output of the 
controller is: 

AP c (t) = -K p .ACE(t) (17) 

Where AP c (t) is output signal of the controller; e is 
constant indicating the specified limit of error signal; 
K P is proportional controller. 

For IACE(t)l<e 

AP c (t) = -K i \ ACE(t)dt (18) 

Based upon the above mentioned facts, the dual 
concept is introduced here in the following way. The 
proportional controller will act during the transient 
period when the error (ACE) is sufficiently larger, 
whereas the integral controller would be the better 
option when the error is small. 

The proposed control scheme is shown in Fig. 6 for 
which, the control law is taken as follows: 



ACE(t) 



► |ACE(t)| 




if |ACE(t)|>£ 




else 



else action 1 



▼ 

-Ki 



if action 

V 

►-Kp 



APc(t) 
— ► 



AP cl (t) = -K pl [ACE l (t)\ 

AP cl (t) = -K a \j ACE x (t)\ 
AP c2 (t) = -K p2 [ACE 2 (t)\ 

AP c2 (t) = -K a \j ACE 2 (t)\ 



for\(ACE l )\>e l (19) 



for\(ACE l )\<s l (2Q) 
for\(ACE 2 )\>s 2 (21) 

for\(ACE 2 )\<e 2 (22) 



Continuous Mode 



FIG. 6 BLOCK DIAGRAM OF THE PROPOSED CONTROLLER 



In the continuous mode, the error signal remains 
within the specified limit, i.e., ACE(t). The integral 
control scheme is the most appropriate for the system 
entering the continuous mode. The frequency and tie 
line power error deviations of hydrothermal system 
under deregulated scenario could be controlled by 
parameters of the dual mode controller (k P , ki and e). 
In both areas, the optimal value of k P , ki have been 
acquired by fuzzy adaptive bacterial foraging method. 

Fuzzy Adaptive Bacterial Foraging Method 

The survival of species in any natural evolutionary 
process depends upon their fitness criteria, which 
relies upon their food searching and motile behaviour. 
The law of evolution supports those species that have 
better food searching ability and either eliminates or 
reshapes those with poor search ability (K. M. Passino, 
2002). The genes of those stronger species get 
propagated in the evolution chain since they possess 
ability to reproduce even better species in future 
generations. Therefore, a clear understanding and 
modeling of foraging behavior in any of the 
evolutionary species, leads to its suitable application 
in any non-linear system optimization algorithm. The 
foraging strategy of E.oli bacteria present in the 
human intestine can be explained by four processes 
namely Chemotaxis, Swarming, Reproduction, and 
Elimination & Dispersal. 

Chemotaxis: The characteristics of movement of 
bacteria in search of food can be defined in two ways, 
i.e. swimming and tumbling together known as 
chemotaxis. A bacterium is said to be swimming if it 
moves in a predefined direction, and tumbling if 
moving in an altogether different direction. 
Mathematically, tumble of any bacteriumcan be 
represented by a unit length of random direction O(i) 
multiplied by a step length of that bacterium C(i). In 
case of swimming, this random length is predefined. 

Swarming: For the bacteria to reach at the richest food 
location (i.e. for the algorithm to converge at the 
solution point), it is desired that the optimum 
bacterium till a point of time in the search period 
should try to attract other bacteria so that together 
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they converge at the solution point more rapidly. To 
achieve this, a penalty function based upon the 
relative distances of each bacterium from the fittest 
bacterium till that search duration, is added to the 
original cost function. Finally, when all the bacteria 
have merged into the solution point, this penalty 
function becomes zero. The effect of swarming is to 
make the bacteria congregate into groups and move as 
concentric patterns with high bacterial density. 

Reproduction: The original set of bacteria, after getting 
evolved through several chemotactic stages reach the 
reproduction stage. Here, the best set of bacteria 
(chosen out of all the chemotactic stages) gets divided 
into two groups. The healthier half replaces the other 
half of bacteria, which gets eliminated, owing to their 
poorer foraging abilities. This makes the population 
of bacteria constant in the evolution process. The 
survival and elimination behavior of any bacterium is 
better known as its motile behavior. 

Elimination and dispersal: In the evolution process, a 
sudden unforeseen event can occur, which may 
drastically alter the smooth process of evolution and 
cause the elimination of the set of bacteria and/or 
disperse them to a new environment. Most ironically, 
instead of disturbing the usual chemotactic growth of 
the set of bacteria, this unknown event may place a 
newer set of bacteria nearer to the food location. From 
a broad perspective, elimination and dispersal are 
parts of the population level long distance motile 
behavior. In its application to optimization, it helps in 
the reduction of the behavior of stagnation (i.e. being 
trapped in a premature solution point or local optima) 
often seen in such parallel search algorithms. The steps 
of the proposed FABF algorithm are described as 
follows: 

Step 1. Initialization 

a) Number of parameters (p) to be optimized. 

b) Number of bacteria (S) to be used to search the 

total region. 

c) Swimming length Ns after which tumbling of 

bacteria will be undertaken in a chemotactic 
loop. 

d) Nc, number of iterations to be undertaken in a 

chemotactic loop(Nc> Ns). 

e) Nre, the maximum number of reproduction to 

be undertaken. 

f) Ned, the maximum number of elimination and 

dispersal events to be imposed over bacteria. 



g) Ped, the probability with which the elimination 

and dispersal will continue. 

h) The position of each bacterium 

P(j,k,l) = [0'j,k,l) I i = 1,2,..., S) at the jth 
chemotactic step, kth reproduction step, and 
1th elimination-dispersal event. 

i) The value of C(i) > 0, where i=l, 2, . . ., S 

denotes a basic chemotactic step size. 

j) The Values Of dattract, Wattract, hrepellant, and Wrepellant 

used for swarming. 

There is a scope to fuzzify the variable C(i) arriving at 
the optimum value of the step size for the given 
problem in less time. Initially, the run length vector 
C(i) value is selected randomly and plays an 
important role in the convergence of SBF 
(Ch.Venkaiah, D.M. Vinod Kumar, 2011) algorithm. A 
small value of C(i) causes slow convergence, whereas 
a large value may fail to locate the minima by 
swimming through them without stopping. The 
selection of C(i) is tedious and time consuming in SBF. 
Hence, fuzzy adaptive scheme is utilized in C(i) to 
ensure the convergence of SBF algorithm. Here, the 
fuzzy input variables are taken as C(i) and the error 
from the objective function to obtain the fuzzy output 
as AC(i) for optimal value. The fuzzy logic rules are 
created by using fuzzy logic tool-box of MATLAB. 

Step 2. Iterative algorithm for optimization: the 
algorithm that models bacterial population chemotaxis, 
swarming, reproduction, and elimination & dispersal 
is given here (initially, j = k = 1 = 0). 

i. Elimination-dispersal loop: 1=1+1. 

ii. Reproduction loop: k=k+l. 

iii. Chemotaxis loop: j= j+1. 

a. For i = 1,2, . . ., S, a chemotactic step for 
bacterium i is taken as follows: 

b. Compute j error aj,k,i) Let 

J error (i,j,k,l) = J„ ror (i,j,k,l)+J cc (6'(j,k,l),P(j,kJ)) 
Let j lmt = j error (i,j,k,l) to save this value since a 
better cost can be drived via a run. 

c. Tumble: generate a random vector A(i) £ R p 
with each element Am(i), m= 1, 2, p, a random 
number on [-1, 1]. 

d. Move: let 

9 i (j + l,k,l) = 6'(j,k,l) + c(i)A(i)p (i)A(O 

This results in a step of size C(i) in the direction of 
the tumble for bacterium i 

e. Compute j ermr (,-, j + i,k,l), and then let 
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j„ ror a, j +u,i) = j emr a, j + ikj) + j cc {o'u + i,k,D, p(j + u,/» 

f. Swim: Let m= (counter for swim length) while 
m< Ns (if it hasn't climbed down too long) 

g. Let m=m+l 

h - If J e m ,r (i, j + h k ,l)< J last (if doing better), let 
lias, = J error d,j + hk,l) and let 

0'(j + \,k,l) = 0'(j,k,l) + c(!)A(i))A r (i)A(i) and use 

this 0'(j + l,k,l) to compute the new 
J error (i,j + l,k,l) as we did in f. Else, let m = Ns. 
This is the end of while statement j. Go to next 
bacterium (i +1) if i * S (i.e. go to b) to process the 
next bacterium. 

i) If j < Nc, go to chemotaxis loop (3) of Step 2. In 
this case, chemotaxis continues, since the life of 
the bacteria is not over. 

ii) Reproduction: J'heaith ='^^ 1 J error (i,j,k,l) 

a. For the given k and 1, and for each i = 1, 2, S, 

J health = V J error (i, j, k, I) be the health of 

bacterium i (a measure of how many nutrients 
it got over its life time and how successful it was 



Initialize all variables. Set all loop 
counters and bacterium index i 
equal to O 



Increase elimination dispersion 
loop counter 
I-I+l 




Increase Reproduction loop 
counter 
k = k+ 1 



at avoiding noxious substances). Sort bacteria and 
chemotactic parameters C(i) in order of ascending 
cost Jheaith (higher cost means lower health) 

b. The Sr bacteria with the highest Jheaith values 
die and the other Sr bacteria with the best values 
split (and the copies made are placed at the same 
location as their parent). 

iii) Chemotaxis loop: j=j+l. If k<Nre, go to 
reproduction loop (2) of step 2. In this case, we 
have not reached the number of specified 
reproduction steps, so we start the next 
generation in the chemotactic loop. 

iv) Elimination-dispersal: for i=l, 2, S, with 
probability p e d, each bacterium is eliminated and 
dispersed (this keeps the number of bacteria in 
the population constant). To do this, if you 
eliminate a bacterium, simply disperse one to a 
random location on the optimization domain. 

v) If l<Ned, then go to elimination-dispersal loop (1) 
of step 2; otherwise end. 

The flow chart of the proposed FABF algorithm is 
shown in Fig. 7. 



Perform Reproduction by killing 
the worse half of the population 
with higher cumulative health and 
splitting the better half into two 




Increase Chemotactic 
loop counter 

j-j+i 



Perform Elimination dispersal 
(For i=l,2,...S with probability 
Ped, eliminate and disperse one 
(to a random location 



Increase bacteriun 
index i=i+l 




Compute the objective function value for 
the i th bacterium as Jerror(iJ,k,l) adding 
the cell to cell attractant effect to nutrient 
(concentration and set Jiast= Jerror(i,j,k,l) 



Tumble (Let the i th bacterium 
take a step of heignt C(i) along 
a randomly generated tumble 
vector A(I» Fuzzify C(i) 



Compute the objective 
function value Jerror(i,j+l,k,l) 
taking into account the cell- 
to-cell attractant effect 




Set Jlast= Jerror(iJ+l,k,l) 

Swim (Let the i th bacterium 
take a step of height C(i) 
along the direction of the 
same tumble vector A(i» 



FIG. 7 FLOW CHART OF THE PROPOSED FABFM 
TABLE 1 COMPARATIVE ANALYSIS 



Mode 


Thermal area 


Hydro Area 


Peak Time 


Overshoot 


Settling Time 


Peak Time 


Overshoot 


Settling Time 


Without SSSC-TCPS-Integral Controller 


3.231 


0.7563213 


6.461 


4.658 


0.9465681 


5.434 


With SSSC-TCPS-Dual Mode 


2.216 


0.0963215 


1.886 


2.365 


0.0975945 


2.864 


% Improvement 


31.41 


87.26 


70.809 


49.227 


89.68 


47.29 1 
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Simulation Result 

A load change of 0.04 p.u M.W in each area has been 
considered to study the comparison between both the 
system with SSSC, TCPS and FABFM based Dual 
Mode controller and the system with conventional 
integral controller and without SSSC, TCPS. A value of 
0.5 has been considered as the gain of integral 
controller. Due to application of FABFM to the system, 
the optimum values of integral controller and 
proportional controllers in both areas i.e., Kii= 
Ki2=7373.26 and K P i=K p2 =8826.125 are obtained. A 
value of ei= £2=0.02 has been considered in this work 
as the specified limit of error signal in the design of 
dual mode controller. Table 1 shows the comparison 
of both the systems. The Discos contract with the 
Gencos as per the following Disco participation matrix: 



DPM 



0.1 

0.0 
0.3 
0.2 
0.2 
0.2 



0.0 
0.1 
0.4 
0.0 
0.3 
0.2 



0.3 
0.0 
0.1 
0.2 
0.0 
0.4 



0.4 
0.2 
0.0 
0.1 
0.1 
0.2 



Table 2 shows the comparison of performance index 
of the system with SSSC-TCPS-Dual Mode and 
without SSSC-TCPS-Integral controller where the first 
one has less performance index than the last one. 

TABLE 2 COMPARISON OF PERFORMANCE INDEX 





Performance Index Value 


With SSSC-TCPS-Dual Mode 


1.086xl0e-5 


Without SSSC-TCPS-Integral 
Controller 


4.7026xl0e-5 



The various frequency deviations and tie line power 
deviations in both areas during a load change of 
0.04p.u MW are shown in Fig. 8. The generations of 
various Gencos in area 1 and 2 are depicted in Fig. 9 
and Fig. 10, respectively. As observed, the system with 
SSSC-TCPS-Dual Mode is far superior to the system 
without SSSC-TCPS-Integral controller in terms of 
peak time, overshoot and settling time in both the 
areas. Three generators in each area have been 
considered for the study. Each Genco participates in 
AGC is defined by following area participation factors 
(apfs): apfl=0.5, apf2=0.25, apf3 =0.25, apf4 =0.5, apf5 
=0.25, apf6 =0.25. Coefficients that distribute ACE to 
several Gencos are termed as "ACE participation 

factors" (apfs). It should be noted that m in ^apfj =1 



.H 



—With SSSC-TCPS-Dual Mode 
— Integral Controller 




30 40 50 
Time [sec] 




— Willi SSSC-TCPS-Dual Mode 
— Integral Controller 



30 40 50 
Timo |sec] 



—Willi SSSC-TCPS-Dual Mode 
— Integral Controller 



30 40 SO 60 70 80 

Time [sec] 

FIG. 8 FREQUENCY AND TIE-LINE POWER DEVIATIONS 
DURING A LOAD CHANGE 



—Willi SSSC-TCPS-Dual Mode 
— Integral Controller 




30 40 SO 

Time [sec] 



-Willi SSSC-TCPS-Dual Mode 
Integral Controller 




pi 

O 

.1 1 

si 

| 0.5 



— Integral Controller 
—With SSSC-TCPS-Dual Mode 




is the number of Gencos. 
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FIG. 10 GENERATION OF GENCOS IN AREA 2 



The comparison between the performance index of 
both the systems is demonstrated in Fig. 11. As seen, 
the system with SSSC-TCPS-Dual Mode has less 
performance index than that without SSSC-TCPS- 
Integral Controller. 

, , — -, 
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FIG. 11 COMPARISON OF PERFORMANCE INDEX 



Conclusions 

In this paper, a methodical FABF technique is 
employed to design dual mode controller for AGC of 
hydrothermal system under a deregulated 
environment. SSSC and TCPS are also used to improve 
the dynamic performance of the system concerning the 
reduction of frequency deviations and tie line power 
deviations during a load change on a two area 
hydrothermal system. This paper focused on the 
development of a new control strategy for SSSC-TCPS- 
dual mode controller is proposed and as well, the 



proposed control algorithm can bring several 
advantages such as follows: 

• The ability to jump out the local optima, the 
convergence precision and speed are 
remarkably enhanced and thus the high 
precision and efficiency are achieved. 

• A control strategy has been proposed to control 
the TCPS phase angle by adjusting the voltage 
of SSSC which in turn controls the inter-area 
tie-line power flow. 

• The first peak frequency deviation of both areas 
and tie-line power oscillations following 
sudden load disturbances in either of the areas 
can be suppressed by controlling the phases 
angle of TCPS or the series voltage of SSSC. 

• The superiority of the proposed SSSC-TCPS- 
dual mode controller is demonstrated owing to 
the less performance index in comparison with 
the system without SSSC-TCPS-Integral 
controller. 

The simulation results reveal that the proposed 
technique successfully mitigates the frequency and tie- 
line power deviations during a load change. 
Furthermore, the tie-line power flow controlled by an 
SSSC units is found to be more efficient and effective 
to improve the dynamic performance of load 
frequency control over inter connected power system 
than that of the system with TCPS or without FACTS 
controllers. The tie line power flow controlled by an 
SSSC can be expected to be utilized as new ancillary 
service for stabilization of frequencies and tie-line 
power oscillations in the congestion management 
environment of the power system. 

Appendix 1 : 

A Fi : The deviation of i th frequency (i=l, 2). 

A P12 : The deviation of the power exchanged between 
the first and second area. 

A PMi : The mechanical output power deviation of 
turbine of i th area. 

A PL : Changing in the electrical load. 

P Ri : The nominal capacity of generation of i* area. 

KRj : The integral coefficient of system of i th area. 

ACEj : The control line of i th area. 

B i : The bias coefficient of frequency of i th area. 
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T Gi : The governor time constant of i th area. 

Ri : The integrating feature of speed (Drop) of area i th 
area. 

T Pi : The time constant of the power system of i th area. 

K Pi : The power coefficient of the i th unit. 

T T : The time constant related to steam inertia in input 
piping and the steam treasury of thermal turbine. 

T RH : The time constant related to steam inertia in the 
process of re-heating of thermal unit turbine. 

F HP : The coefficient related to power generated in the 
high pressure steam part of thermal turbine. 

F LP : The coefficient related to the power generated in 
the medium pressure steam part of thermal turbine. 

T w :The time constant of initiating water in hydro 
turbines. 

R F : The integrating speed feature of continual state of 
mechanical governor in each unit. 

R T : The integrating speed feature of transient state of 
mechanical governor in each unit. 

T R : The reset time constant of mechanical governor in 
each unit. 

Kcp: The coefficient of phase shifter system. 
T F : The time constant of phase shifter system. 
/ : System objective function. 

Appendix 2: 

• Data for Thermal Reheat Power System (LA. 
Chidambaram and S. Velusami, 2005) Rating of 
each area = 2000 MW, Base power = 2000 MVA, 
fo = 60 Hz, Rl = R2 = R3 = R4 = 2.4 Hz / p.u.MW, 
Tgl= Tg2 = Tg3 = Tg4 = 0.08 sec, Trl = Tr2 = Trl 
= Tr2 = 10 sec, Ttl = Tt2 = Tt3 = Tt4 = 0.3 sec, 
Kpl = Kp2 = 120Hz/p.u. MW, Tpl = Tp2 = 20 
sec, fi x = fi 2 = 0.425 p.u. MW/Hz, Krl = Kr2 
=Kr3 = Kr4 = 0.5, In T12 = 0.545 p.u.MW / Hz, 
al2 = -l, A roi =0.01 p.u.MW 

• Data for TCPS (R. J. Abraham, D. Das and A. 
Patra, 2006): K<j> =1.5 rad/Hz, ^ max = 10° , 
0min = -10", Tps = 0.1 sec. 

• Data for SSSC (I. Ngamroo, 2001): Tsssc=0.03 
sec, Tw=10 sec 
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